Coherent quantum effects through dispersive bosonic media 
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The coherent evolution of two atomic qubits mediated by a set of bosonic field modes is investi- 
gated. By assuming a specific encoding of the quantum states in the internal levels of the two atoms 
we show that entangling quantum gates can be realised, with high fidelity, even when a large number 
of mediating modes is involved. The effect of losses and imperfections on the gates' operation is 
also considered in detail. 
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I. INTRODUCTION 

There is great emphasis, in current research, on the 
identification of physical systems and conditions where 
coherent quantum effects are susceptible to emerge. In 
fact, the achievements of quantum communication and 
key distribution quantum metrology 0, and the 
prospects opened by quantum computation [J, Q are all 
centred on quantum coherence. Also, the latter has been 
shown to occur in biological proteins and is suspected 
to play a role in some mechanisms of energy transfer in 
biological structures Q. In this context, it is extremely 
relevant, for both fundamental and applied reasons, to 
single out systems and operating regimes where coherent 
effects can take place, and to analyze them in detail. 

A recurrent paradigm for quantum information pro- 
cessing comprises a set of static qubits where the in- 
formation is stored (usually embodied by internal levels 
of atoms, ions or quantum dots in solid state devices), 
interacting through a bosonic field (light or vibrational 
phonons). Countless variations of this paradigm have 
been proposed and studied, in both theory and practice, 
over the last fifteen years, with peculiar advantages and 
disadvantages for the realization of controlled coherent 
manipulations fol— T30l| . 

However, in all of these studies, only one of the bosonic 
field modes - or a selected few, in more sophisticated set- 
tings - are assumed to mediate the interaction between 
the qubits. This assumption may be reasonable if the 
space where the bosonic field modes resonate - which we 
will henceforth refer to as the "fibre" - is small, so that 
the frequency spacing is very large, and if the experi- 
mental control is very accurate. Such conditions are not 
always met in practical situations: often more than one 
mode mediates the interaction between the qubits, for in- 
stance when the fibre is too long, or when unwanted de- 
grees of freedom outside direct experimental control enter 
the dynamics. When such situations occur, one expects 
the capability of the system to mediate the controlled, co- 
herent evolutions for the two qubits to be degraded, sim- 
ply because the quantum information disperses among a 
greater number of degrees of freedom. The question then 
arises of whether controlled dynamics arc still possible in 
"dispersive" bosonic media. 



This paper intends to launch an investigation to un- 
derstand whether, to what extent, and under what con- 
ditions may coherent quantum effects be mediated by 
multiple degrees of freedom. Notice that such a ques- 
tion is relevant not only to quantum information pro- 
cessing (where, it may be argued, the relative flexibility 
in the engineering of quantum systems sometimes allows 
one to circumvent the effect of undesired degrees of free- 
dom), but also to other, more fundamental applications 
of quantum theory, like in the modeling of certain light- 
harvesting biological compounds, where several frequen- 
cies of the light field are thought to coherently mediate 
the transfers of energy. 

In this study, we shall refer to a specific basic 'proto- 
type' which has drawn considerable attention as a can- 
didate for 'distributed' quantum information processing. 
Two distant atoms arc trapped in two remote optical 
micro-cavities, see Fig.Q] and interact locally with one of 
the cavity modes. The cavities are then in turn connected 
through an optical fibre where many light modes couple 
to the cavity modes. In practice, this prototype might 
represent two atoms or ions trapped in micro-fabricated 
optical cavities directly coupled to the same integrated 
optical fibre [3lj . However, our results will apply to more 
general systems sharing the same formal description. In 
this respect notice that, by transforming to normal modes 
of the cavities plus fibre field, our system can be shown 
to be equivalent to two qubits whose interaction is di- 
rectly mediated by a set of bosonic modes (without the 
two interfacing cavities fields). 

To demonstrate that controlled coherent quantum evo- 
lutions arc possible, we shall focus on the realisation of a 
quantum controlled-z (CZ) gate for the two qubits. The 
perfect implementation of such an operation would allow 
for maximal entanglement to be generated and, supple- 
mented with single-qubit gates, would provide one with a 
universal set of quantum gates, on any number of qubits 



FIG. 1. (color online). The prototype: two remote nodes, 
where atoms are trapped in, are linked by a media containing 
many bosonic modes at different frequencies. 
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system. On the other hand, wc want to identify a situa- 
tion where a controlled entangling quantum gate can be 
mediated by several degrees of freedom, so that we have 
to find a way to reduce the effects of dispersion through 
such a medium. 



Both these issues will be solved by encoding the two 
qubits asymmetrically - i.e., differently in the two sites 

in one excited state and two kinds of ground states, 
following a strategy introduced in Ref. [29|]. In the first 
cavity, labelled by 1, the qubit is encoded in the ground 
state | g)\ and in the excited state [ e)i, which are coupled 
to each other via the cavity mode, by a rotating-wave 
Hamiltonian. In the second cavity, labelled by 2, the 
ground state \g) 2 is also coupled to the excited state 
| e) 2 , but the qubit is encoded in the ground states | g) 2 
and I s) 2 , and the latter is not coupled to any state and 
does not evolve at all in the dynamics we consider. This 
configuration allows us to study the realisation of the CZ 
gate by restricting to the single excitation subspace, thus 
both reducing the computational effort required — if iV is 
the number of mediating modes, the size of the relevant 
Hilbert subspace only scales like N - and limiting the 
dispersive effect of the modes. 

We shall assume the cavity fields ai and 02 to be at res- 
onance with the atomic transitions, so that the rotating- 
wave approximation holds, and shall have them interact- 
ing with N modes of the "fibre" at different frequencies. 
In the frame rotating at the cavity frequency, the total 
Hamiltonian of the system is then 



N I 2 \ /k=2.j=N \ 

H = Aw J b % + 9 kQk I l fe +h - c - + £ Vkb i a l + h - c - ' w 
3=1 \fe=i / yfe=i,3=i J 

where Awj is the frequency difference between the kth mediating mode and the cavity mode, a\ (a*) and bj(bj) are 
the creation (annihilation) operators for the cavity modes and the fibre modes, respectively, v k denotes the coupling 
strength between the mode of cavity k and the mediating modes (possibly including a phase), and g k represents the 
coupling strength between atom and field in cavity k. Note that the fibre-cavities couplings depend on the cavity 
but are assumed to be the same for all the fibre modes. The emphasis in our analysis is rather on the effect of the 
detunings A.Wj of the fibre modes. 

Also, we will consider both losses of the cavities and fibre and spontaneous emission from the atoms (by far the 
main sources of decoherence at optical frequencies), so that our dynamics is governed by the master equation 

2 2 TV 

p = -i[H, p) + k ]T L[a k ]p + 7 ^i[ ( T-]p + r^ L[b k ]p , (2) 
fe=i fc=i 3=1 

where L[d] = 2opo* — 6^ dp — pd^d is for operator 6, and 7, k and T are the atomic spontaneous emission rate and the 
loss rates of the cavities - assumed, for simplicity, to be identical - and of the fibre, respectively. 

Initially, all the field modes including the cavity and mediating modes are in the vacuum state | 0)/. Moreover, to 
check the reliability of the controlled phase gate, we consider the generic input state of two qubits 

= a I e)i I g) 2 + f3 \ e) 1 \ s) 2 + 7 | g)i \ g) 2 + S \ .9)1 | s) 2 , (3) 

where a, /3, 7 and S are distributed according to the Haar measure of U(A) (that is, as resulting from the application 
of a random Haar-distributed unitary on any fixed normalised state). Let us recall that the Haar measure is defined as 
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FIG. 2. (color online). The atomic level configuration of the 
two atoms (k — 1, 2). Notice that level | s)i does not play any 
role in our scheme. Level \ e)2, instead, enters the dynamics 
but is not used to encode any quantum information. 



[32j. No specific initial state will be assumed: the fi- 
delity of the gate will be averaged over the unitarily- 
invariant Haar measure of the two-qubit state space. By 
adopting a particular structure for the atomic levels and 
qubit encoding, we will show that, in several situations, 
this non-trivial controlled quantum dynamics of the two 
qubits can be approximated with high fidelity, even if 
more than one mode mediates the interaction between 
the distant systems. This also implies that the creation 
of nearly maximal entanglement between the two qubits 
is possible. 



II. METHODS 

The main obstacles to our programme are twofold. On 
the one hand, we have to deal with the customary com- 
putational difficulty of describing a many-body quantum 
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the measure which is invariant under both left and right multiplication by any unitary transformation. A distribution 
borrowed from the Haar measure is a natural choice for pure quantum states if one does not want to privilege any 
specific direction in the Hilbert space, and hence for testing the average fidelity of a given quantum operation, as is 
the case here. The desired, ideal output state | ^>) ut of the CZ gate, corresponding to the input state | ^)i n , is 

V)out = a | e)i | g) 2 - p | e)i | s) 2 + 7 I g)i I 9)2 + 6 \ 9)1 | s)a , (4) 
where the phase flips only for the state | e)i | s) 2 - 



Clearly, of the four superposed states defining | ^)i n , only 
I e) 1 1 g) 2 and | e) 1 1 s) 2 evolve, which further simplifies our 
task (besides the fact that only the single excitation sub- 
space is involved) . It is therefore straightforward to inte- 
grate Eq. @ obtaining the final state g(t) of the system 
given the initial state | ^)i n \ 0) /, and to evaluate the Haar 
average of the fidelity: 

F = [ U T17 [g(t)\ I tf ) out d| *} m , (5) 

J Haar 

where the notation J Haar d| ^)i n loosely refers to the fact 
that integration is carried over Haar-distributed input 
states, and Tif stands for the trace over the field's de- 
grees of freedom. In practice, this integral has been car- 
ried out by sampling the input states according to the 
Haar distribution and by averaging the resulting fidelity. 
Samples of 200 points turned out to give reliable esti- 
mates, in that the average fidelity at the first two signifi- 
cant figures would not change by increasing the sample's 
size. Henceforth, we will always report fidelities at two 
decimal figures. For some finer comparisons, we improved 
our sensitivity by increasing the integration samples. In 
these cases as well we will limit ourselves to two decimal 
digits, and will just qualitatively point out the configu- 
rations providing higher fidelities. 

III. RESULTS: CONTROLLED COHERENT 
EVOLUTIONS 

We will study now the approximated realization of the 
CZ gate as detailed in the previous section, in several 



different situations, mainly varying the number N of fi- 
bre modes and their detuning with respect to the cavity 
modes. All the results obtained are summarized at the 
end of the section. 

Hereafter, we refer to the fibre mode resonant with the 
atomic transition between | e)k and | g)k as the 'central 
resonant mode'. In order to find regimes with high av- 
erage gate fidelity, we set g\ = V\ = g, while changing 
52, v 2 and the interaction time t. To make comparisons 
more clear, we define two parameters Sg = (<?2 — g) / g 
and 5v = (v 2 - g)/g- 



A. Two mediating modes. The Minimum Gap. 

To study coherent effects resulting from the competi- 
tion of multiple mediating modes, we start from the case 
of only two mediating modes, with the same absolute 
detuning A with respect to the frequency of the central 
resonant mode. Fig. [3] shows the average gate fidelities 
F for different detunings at a time where the first peak 
in fidelity is achieved. Clearly, F decreases as the detun- 
ing A increases, since the two fibre modes become more 
and more off-resonant. For instance, at the peaks, F is 
0.99, 0.98 and 0.96 for A/g = 0.1,0.2,0.3, respectively. 
To put these fidelities into some context, let us mention 
that, depending on the noise and on the affordable com- 
putational overhead, error thresholds between 10 -4 and 
10 -2 are believed to be required to achieve fault-tolerant 
quantum computation (33j . 

Typically, when A reaches 0.45g, and the gap between 



FIG. 4. (color online). The average gate fidelity F for two 
detuned modes versus Sg and Sv for loss rates 7 = n — V = 
W~ 2 g, and for A = O.lg and gt = 4.6 (a), A = 0.2g and 
gt = 4.56 (b), and A = 0.3g and gt = 4.54 (c). 

the two modes becomes comparable to the interaction 
strength g, the average gate fidelity F drops to 0.90. As 
a reference, we will say that 0.9<? is the "minimum gap", 
within which "high fidelity" (90%) can still be recovered 
with two mediating off-resonant modes. In a cavity QED 
implementation with g ~ I GHz, this gap would corre- 
spond to a fibre length I of around one meter, since the 
spacing of two neighbouring mediating modes is approx- 
imately cir/l ~ 10 9 Hzm/Z. 

As shown in Fig. [4j F decreases by around 4%, when 
losses are included for 7 = re = T = 10~ 2 g. Further nu- 
merical analysis showed that for 7 = re = V = 10~ 3 g, 
the optimal fidelities decrease by around 0.5%. These 
findings are in agreement with previous studies [2(| ar- 
guing that, for nearly resonant couplings like the present 
one, average fidelities are virtually unaffected for loss to 
coupling strengths ratios around 10 -4 , are only slightly 
affected for ratios around 10~ 3 , and start to affect the 
fidelity significantly - between 1 and 10% - when such 
ratios reach 10~ 2 . Of course, observing nearly ideal lev- 
els of coherence requires a high degree of isolation. How- 
ever, notice that a quantum CZ gate implemented with 
a fidelity of 85% would still be a remarkable signature 
of quantum coherence, implying the creation of substan- 
tial entanglement and of coherent off-diagonal elements 
in the density matrix of the two qubits. 

A note about the degree of stability of the gates ob- 
tained is also in order here: a variation of the order of 
0.1 in the parameters Sg and Sv, corresponding to a vari- 
ation of about 10% in the coupling strengths, leads to 
a decrease in the gate fidelity of around 1%. The same 
degree of sensitivity occurs with respect to fluctuations 
in the interaction times. Notice that, while achieving 
a lower peak in fidelity, the performance of more off- 
resonant modes (A = 0.3g in the figures) is slightly less 
sensitive to imperfections, as evident from Fig. [3] 



B. Increasing number of mediating modes. 

Let us now consider the case of a larger number of 
mediating modes. Fig. [3] (a) shows the case of 31 medi- 
ating modes, with frequencies equally and symmetrically 
spaced around the central resonant mode and filling the 



(a) (b) 

FIG. 5. (color online). The average gate fidelity F for a 
set of equally spaced mediating modes within the minimum 
gap A = 0.45<? versus Sg and Sv at the time 4.3g~ 1 , without 
considering losses. In (a) 31 modes are considered; in (b) the 
central resonant mode has been removed. 

minimum gap. At the peak, a remarkable average fidelity 
F = 0.99 can still be achieved: the competition between 
the distinct modes, which could favour the 'dispersion' 
of the coherence, is clearly canceled out at optimal times 
for such a configurations of closely packed, nearly reso- 
nant modes. In the cavity QED model, for g ~ 1 GHz, 
the spacing considered here (0.03g) would correspond to 
a fibre of approximately 30 m: in principle, very long res- 
onators can still mediate quantum evolutions coherently. 
In Fig. [5] (b) the same case without the central resonant 
mode is depicted: at the peak, the average fidelity drops 
very slightly but is still well above 0.99, which proves 
that the coherent evolution mediated by these 31 modes 
is not an effect due to the presence of the central res- 
onant mode. Nor is this high fidelity a consequence of 
the symmetric distribution of the modes around the cen- 
tral resonant frequency: this has been directly tested by 
shifting all the frequencies of the mediating modes and 
does not produce any significant alteration in the maxi- 
mal average fidelity, as apparent from Fig. [5] 

To evaluate the influence of the modes very close to 
resonance with the cavity frequencies, we now turn to 
cases where the set of mediating modes are all outside the 
'minimum gap'. In Fig|6] (b), 32 equally spaced modes 
are considered: half of them spans the range between 




"(a)" (b) 



FIG. 6. (color online). The average gate fidelity F for equally 
spaced modes outside the minimum gap. with detuning rang- 
ing from 0.45g to 0.9g, and from — 0.45<7 to —0.9g versus Sg, 
at the time 4.48g -1 , without considering losses. In (b), 32 
modes are considered; in (a), a 33 th mode has been added, at 
the central resonant frequency. 



FIG. 7. (color online). The average gate fidelity F for equally 
spaced modes outside the minimum gap, with detuning raging 
from 0.45(7 to 0.9g, and from — 0A5g to —0.9g, at the time 
4.48<? _1 , considering loss rates 7 = k = F = 10~ 2 g. In (b), 32 
modes are considered; in (a), a 33 th mode has been added, at 
the central resonant frequency. 



(a) (t» 

FIG. 8. (color online). The average gate fidelity F for 16 
equally spaced mediating modes outside the minimum gap, 
with detuning ranging from 0.45g to Q.9g, and from —0A5g 
to — 0.9<7 at the time 4.48(? _1 . No losses are considered in (a), 
while loss rates 7 = k = T = W~ 2 g are introduced in (b). 



0.45(7 and 0.9g, while the other half spans the range be- 
tween —0.45(7 and —0.9(7. Even in such an off-resonant 
configuration, a maximum average fidelity of 0.98 can be 
achieved. This is a non-trivial finding, mostly if com- 
pared to the - much lower - fidelity achievable with only 
two modes with frequencies at the minimum gap (which 
is about 0.9). In this case, more modes, and farther 
off from resonance, allow one to achieve a better per- 
formance in terms of mediated quantum coherence. This 
is therefore a remarkable instance where the cooperation 
between the mediating modes prevails over their compe- 
tition, and their effect could hence be harnessed to im- 
plement distributed coherent quantum evolutions. Fig. [6] 
(a) confirms that the addition of the central resonant 
mode has only a relatively modest impact on the max- 
imal fidelity (which raises slightly but is still essentially 
0.98). 

Fig. [7] displays the effect of losses on the optimal av- 
erage fidelity: for decay rates 7 = k = T = 10~ 2 g, the 
average gate fidelity F drops from around 0.98 to around 
0.94. Concerning decoherence, the resilience of many 
modes seems to be comparable to that of few mediat- 
ing modes. 

Finally, in Fig. [5J the case of 16 modes with detun- 
ing ranging from 0.45 g to 0.9g is considered. The results 
are very similar to the cases of 32 and 33 symmetrically 
detuned modes (including center resonant mode), in the 
same regime and interaction time. At the peak, the av- 
erage gate fidelity is 0.99. Such a fidelity is reduced to 
0.95 when losses with 7 = k = T = 10~ 2 g are taken into 
account [see Fig. |5](b)]. 

The situations addressed above, where finite regions 
of frequencies are not populated, although of clear the- 
oretical interest in the context of our study, might seem 
rather artificial in practice. However, even for bosonic 
fields, such situations could be of relevance in systems 
like photonic crystals, where photonic bandgaps arise for 
properly modulated refractive indexes [34j . 



C. 'Many' mediating modes. 

Raising the number of mediating modes to 100, one 
can see from Fig [10] (a) that very high fidelities can be ob- 
tained, still above 0.99, if all the modes are taken within 
the minimum gap, which is promptly explained by the 
presence of a large number of modes very close to res- 
onance. More interestingly, in Fig. [TU] (b) we have re- 
considered the situation of Fig. [5] (b) , by adding to it 35 
modes in each direction with respect to the central reso- 
nant frequency. These additional modes would of course 
be present in the realistic modelling of a 300 m long fi- 
bre, and could approximately account for the whole field 
resonating in the fibre (as the effect of more and more 
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FIG. 9. (color online). The average gate fidelity F versus time 
for Sg = 0.9 and Sv = 0, when the frequency spacing is kept 
constant but the centre of the set of frequencies is shifted. 
In (a), the frequency difference between neighbouring modes 
is 0.2g and the central frequency is shifted from to 0.2g; 
the curves at the highest peak refer to two mediating modes, 
while the two curves at the lowest peak refer to 30 mediating 
modes. In (b), the frequency difference between neighbouring 
modes is 0.9g and the central frequency is shifted from to 
0.9p; the curves at the highest peak refer to two modes, while 
the two curves at the lowest peak refer to 30 modes. No- 
tice that, mostly around the peaks, the curves are essentially 
indistinguishable from each other. 



6 




(a) " (b) 



FIG. 10. (color online). The average gate fidelity F for a set 
of 100 equally spaced mediating modes versus 5g and Sv at 
the time 4.3<? _1 , without considering losses. In (a), the modes 
are within the minimum gap, with detuning A from — 0.45<? 
to 0.45<?. In (b) the modes have detuning A ranging from 
— 1.5<? to 1.5g. The plot (b) depicts the same case of Fig.[S|b) 
but with the addition of 70 modes outside the minimum gap. 

off-rcsonant modes can be considered to be very small). 
In this case as well, the optimal average fidelity has been 
found to be well above 99%. 

As for the stability of the quantum operations in the 
face of imperfections, a comparative analysis of the plots 
shows that a larger number of modes grant flatter and 
flatter fidelity peaks [compare, in particular, the two 
peaks of Fig. EJa-b)]. So, while the fidelity obtained is 
slightly lower, the stability allowed by many mediating 
modes is higher than that of two, or even one (20| . medi- 
ating modes (it can be inferred from Fig.[T0]that, for 100 
modes, a 10% variation in the coupling strengths leaves 
the fidelity practically unchanged at two decimal figures). 

D. Summary of results. 

Summing up, in this study we have shown that: 

• for only two mediating fibre modes, the optimal av- 
erage fidelity F of a CZ decreases relatively rapidly 
with increasing gap between the modes, plummet- 
ing to 0.9 for a gap equal to 0.9<? (where g is a 
reference value for the atom-cavity and cavity-fibre 
couplings); 

• increasing the number TV of mediating modes, one 
can still achieve highly reliable gates. For TV ~ 
100, F ~ 0.99. This is the case regardless of the 
symmetry of the frequency spacing of the mediating 
modes with respect to the central resonant mode, 
and of the presence of a mode at resonance with 
the cavities; 

• nearly perfect gates (F ~ 0.98 for N ~ 30) can even 
be obtained for off-resonant interacting modes. In 
this large number of modes actually outper- 
forms fewer modes; 

• loss rates around 10~ 3 g are necessary to operate 
such gates in perfect conditions (to all practical 
purposes); however, even rates around 10~ 2 g allow 



for coherent effects to emerge (and entanglement to 
be generated); 



• for two mediating modes, imperfections of the or- 
der of 10% in the coupling strengths or interaction 
times affect the resulting fidelity by approximately 
one percent; this stability improves for increasing 
number of modes, and the fidelity is virtually un- 
changed for imperfections around 10% and 100 me- 
diating modes. 



IV. CONCLUSION AND OUTLOOK 

By adopting a specific qubit encoding, and meaning- 
fully restricting the system to the single excitation sub- 
space, we gave stringent evidence that distributed co- 
herent evolutions can be mediated by a large number of 
bosonic modes. 

In practice this means that, if losses can be kept at 
bay, 'very long' resonating structures could be adopted 
to mediate distributed quantum computation (or, con- 
versely, that relatively low coupling strength could still 
be admissible). In some particular configurations of the 
mediating frequencies, where a bandgap is present, wc 
also found cases where the collective effect of many modes 
outperforms the mediation of fewer degrees of freedom, 
thus dispelling the notion that employing fewer degrees 
of freedom where, roughly speaking the 'quantum infor- 
mation is less likely to disperse', should result in better 
performances (it should be noted though, to avoid any 
ambiguity, that all the fibre degrees of freedom are actu- 
ally to be considered as 'under control'). 

At a more fundamental level, this study shows that 
coherent quantum effect can be mediated by many modes 
and, potentially, over long distances. 

As a future step, we aim to extend this investigation 
beyond the restricted, specific scheme we adopted here, 
to more general situations where distributed quantum 
evolutions are mediated by a number N bosonic modes. 
We intend to address situations, whose building block is 
essentially the "spin-boson" model [35[ , where N is large 
enough that it cannot be treated with elementary tech- 
niques, but small enough not to allow a reliable descrip- 
tion of the field in terms of master equations for common 
baths, like in (36l . Let us mention that, very re- 

cently, the ground state properties of a "two-impurity" 
spin-boson model have been investigated by adopting a 
variational approach [38|, and t-DMRG methods have 
been adapted to the study of dimers strongly interact- 
ing with bosonic baths [39||. To tackle the study of the 
dynamics of such models, we will instead resort to tech- 
niques based on co-moving bases and Ehrenfest dynamics 
[40j. This investigation is currently under way [4l|. 
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